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LINE BROADENING OF THE TRIPLET DIFFUSE SERIES 
OF HELIUM IN AN ARC-JET PLASMA 


Carl Douglas Scott, B. A. 

The University of Texas at Austin, 1969 

Supervising Professor: Dr. L. W. Frommhold 

Abel-inverted spectral-line profiles of the helium lines in the 

3 3 

2 P-n D series were measured in an arc-jet plasma flow that had an electron 

15 —3 

density of 1.4 x 10 ^ cm and a temperature of 7800° K at the center line. 

The electron density and temperature radial distributions were found from 
Abel-inversions of absolute continuum intensities and relative line inten- 
sities, respectively. The measured profiles were compared with theoretical 
profiles which were calculated from the quasistatic Stark broadening theory 
(similar to the theoretical calculations of Pfennig and Trefftz). In these 
calculations, Doppler broadening, although negligible, was included as a 
device to avoid singularities. An empirical correction to the theory was 
made which permits the validity of the quasistatic calculation for electrons to 
be extended into the impact regime toward the line cores and which improves 
the agreement with the measured profiles. 
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1 . 0 INTRODUCTION 


The study of the spectra of radiation emitted from plasmas has long 
been an invaluable tool in atomic research, not only for determining proper- 
ties of the constituent atoms in the plasma, but also for determining proper- 
ties of the plasma as a whole and for determining the effects of the plasma 
on the atoms. In this work, established spectroscopic techniques were used 
to determine the electron density and temperature in a helium arc-jet plasma. 
Then, these measured parameters were used to calculate the Stark broadening 
of spectral lines in the triplet-diffuse series of helium. These calculations 
were made by use of the quasistatic approximation for ions and electrons, and 
were compared with measured profiles in a plasma where the center-line elec- 
tron density was 1.4 * 10 ±13 percent cm and the center-line temperature 

was 7^00 ± 17 percent °K. Not only is knowledge of the broadening useful 
from a laboratory viewpoint, but it is useful in determining conditions in 
stellar plasmas, and is important in solving radiative-transfer problems. 

The theory of Stark broadening has been worked out for two limiting 
approximations: the quasistatic approximation, which applies to low-density 

and low-temperature plasmas; and the impact approximation, which applies 
to high-density and high-temperature plasmas. 

The application of these approximations to helium line broadening has 
been restricted until now to comparisons with experimental measurements only 
in the extremes where one theory would be expected to predominate over the 
other. There has been increased astrophysical interest in the helium-triplet 
diffuse series because it is strong in many stellar spectra. This interest 
is also in conditions between the extremes, that is, in a range of electron 
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densities of lO 1 ^ to lO 1 ^ cm" 3 . This is the range of conditions obtainable 
in arc- jet flows. 

Measurements of the profiles of helium lines have been predominantly 
made in either the relatively low-density plasmas of radio-frequency (rf) dis- 
charges or the relatively high-density plasmas such as arcs and shock tubes. 

Wulff (ref. l) measured the profiles of various lines in the sharp, 
principal, and diffuse series of neutral helium in a pulsed arc discharge 

where the electron density was 3 x 10 and the temperature was 30 000° K. 

° o 

Comparison of the 4471 A and 4026 A helium lines with a statistical theory of 
the Balmer and H_^ lines, respectively, did not yield good agreement, 

as might be expected, because helium has a much more complex term structure 
than hydrogen does. Griem, et al. (ref. 2) calculated the helium line broad- 
ening in the impact approximation. They compared these calculations with 

o o 

those measured by Wulff (ref. 1) 4713 A and 3965 A lines, and found good 
agreement . 

Vidal (ref. 3) made accurate measurements of the line broadening of 

the singlet and triplet diffuse series of helium in an rf discharge. In 

13 -3 

Vidal's experiments, the electron density was 2.65 x 10 cm , and the tem- 
perature was 1850° K. Pfennig and Trefftz (ref. 4) calculated the line wing 
profiles by using the quasistatic theory and then they compared these pro- 
files with the measurements (ref. 3). Agreement in the far wings was 
reasonably good. 

0 3 3 

Greig, et al. (ref. 5) measured the profiles of the 3889 A (2 S-3 P) 
and the 5016 A (2 S-3 P) helium lines from a high-density plasma in a shock- 
tube source. The electron density determined from the width of the line 

was in the range of 10"*"^ to 10"*"^ cm 3 . They found that the widths of the 



3 


helium lines just discussed compared well with the generalized impact theories 
of Griem, et al. (ref. 2) and Oertel (ref. 6). 

In these experiments, measurements of line profiles in the helium- 
triplet diffuse series were made in the free-stream flow of an arc- jet plasma. 
The electron density was determined from absolute continuum-intensity measure- 
ments , and the electron temperature was determined from relative total line 

intensities. In the present work, an attempt is made to use the quasistatic 

3 3 

theory to describe helium 2 P-n D line profiles and to include an empirical 
correction to extend the range of validity of the theory toward the center of 
the line. This correction is patterned after those who applied it to hydrogen 
broadening, Schluter (ref. 7) and Edmonds, et al. (ref. 8). This correction 
makes the perturber density N a function of the wave number Av, tending 
from N = N g in the core to N = 2N^ in the far wings. 

Because of the inhomogeneous nature of the free-stream jet, the meas- 
ured profiles were transformed to profiles which are functions of the jet 
radius by Abel-inversion. Comparisons of the experimental results and the 
theory are made at the center-line radius r = 0. The temperature and the 
electron density were determined as functions of the plasma radius by Abel- 
inversion of line and continuum intensities, respectively. The theoretical 
bases for these calculations are described in Section 2.1, Line Intensities 
and Temperature, and Section 2.2, Continuum Intensity and Electron Density. 

The theoretical line-profile calculation in which the Doppler profile 
is folded with the quasistatic profile is described in Section 2.3, Quasistat- 
ic Stark Broadening Theory. This folding primarily is a mathematical conven- 
ience because the Doppler broadening is negligible compared with the Stark 
broadening at these conditions. The quasistatic profile is calculated from 
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the relative line strengths and frequency shifts of the Stark components of 
Pfennig and Trefftz (ref. 9). Baranger and Mozer's high-frequency component 
field-strength-distributi on function, applicable for the electrons, was used. 
The ions are also accounted for by the same field-strength-distribution func- 
tion. The empirical correction is also described. 

The facility and instrumentation are described in Section 3.1, The 
Apparatus; the spectroscopic measurement techniques are described in Sec- 
tion 3.2, The Spectroscopic Measurements. These measurements include the 
Abel-inverted line and continuum intensities, Abel-inverted line profiles, 
and the instrument profiles. 

The results of the measurements and the theoretical calculations are 
discussed in Section 4.0, Discussion of Results and Comparison of Experiment 
With Theory. It is shown that the agreement in the wings is reasonably good, 
except in regions where the forbidden components contribute to the profiles. 
The empirical correction is shown to improve agreement toward the line cores 
in cases where errors in the measured profiles are not too great. 

The summary is given in Section 5-0, Summary. The relations involved 
for relative and absolute-line and continuum-intensity calibration are shown 
in appendix A. The cases for broad and narrow lines are discussed. The Abel- 
inversion calculations are discussed in appendix B. 



2 . 0 THEORY 


2.1 LINE INTENSITIES AND TEMPERATURE 

Temperature is a quantity which defines the distribution of energy in 
a statistical system in equilibrium; hence, one can determine the temperature 
i f the population of energy levels in an electronic system is known. Because 
the intensity of line radiation emitted from an optically thin gas is propor- 
tional to the population density of an energy state, one can obtain the tem- 
perature from intensity measurements. 

The electron temperature in the arc-jet plasma was determined by the 
use of the standard method of the Boltzmann plot, wherein relative total-line 
intensities are used and log (iA^/gf) is compared with E^; I is the rela- 
tive total intensity, X is the wavelength of the line, g is the statistical 
weight of the lower energy level of the transition, and f is the absorption- 
oscillator strength. The slope of the resulting line equals -(l/kT), where 
k is the Boltzmann constant and T is the electron temperature. This method 
is based on the assumption that the involved lines are optically thin. Also, 
the method is based on the expression for the intensity of a spectral line 

n 

2irhc r -E /kT 

I = 3— ^ gfN Q e U/ (1) 

X 

where h is Planck's constant, c is the velocity of light, r Q is the 
electron radius, and N q is the total atom density. This equation is based 
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upon a local thermodynamic equilibrium, because a Boltzmann distribution of 
the energy levels is assumed. 

The assumption of optical thinness is implicit here because the radi- 
ance N(r) is determined from an Abel-inversion which holds for an optically 
thin plasma. The Abel-inversion numerical technique is discussed in appendix B. 

2.2 CONTINUUM INTENSITY AND ELECTRON DENSITY 

To determine the electron density, the equation for the continuum in- 
tensity was solved for the electron density. The equation used was that given 
by Anderson and Griem (ref. 11 ) for the emission coefficient for the continuum 



where a is the fine structure constant, a is the Bohr radius , c is the 

o 

speed of light, E is the ionization energy of hydrogen, k is the Boltzmann 
constant, T is the temperature, A is the wavelength, n is the principle 
quantum number , N g = is the electron density, and N^/N^ is the rauio of 
ion ground-state population to total ion population (set equal to l). 

The relative photoionization absorption cross section for the n,£ 
levels is g , in units of the Kramers-Unsold cross sections, but is multi- 
plied by the ratio of the statistical weight of the n,£ level to that of the 

2 

hydrogen levels 2n and is divided by the statistical weight of the ion 
ground state. The variable g n (A) is the bound-free Gaunt factor, and g (A,T) 
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is the free- free Gaunt factor averaged over a Maxwellian velocity distribution 
of electrons. The values of these Gaunt factors were interpolated from the 
tables of Karzas and Latter (ref. 12). The E m is the ionization energy; E ^ 
is the excitation energy of the level. 

The variable AE is the advance of the series limit given by Griem 
(ref. 13) in the formula 


AE = 


1 + 



(3) 


The variable AE^ is the reduction in ionization energy given by Griem 
(ref. 13) in the formula 


AE 


00 




(h) 


where E is the photon energy and e is the permittivity of free space. 

In the sums , only those terms contribute for which the photon energy 
is greater than the energy required to excite a bound electron into a free 
state, taking into account the advance of the series limit. The upper limit 
on n for the bound- free sum is one less than the n corresponding to the 
last discernible line in the series. The electron density is then computed 
from the absolute continuum-intensity measurement 


N (r) 
e 


N(r,A) 


:(X) 


1/2 


( 5 ) 


where N(r,A) is the power radiated per unit volume per steradian per unit 
wavelength at a point in the plasma at a distance r from the center line. 



2.3 QUASISTATIC STARK BROADENING THEORY 


The theory of Stark broadening of spectral lines may be approached 
from two opposing idealizations: the impact theory and the quasistatic theory. 

The impact theory is applicable when the phase of an emitted wave is 
on the average, only slightly changed during the time of the order of a few 
collision times. This condition applies in two general cases. In the first 
case, the collisions are strong enough to disrupt the phase of the wave, are 
well separated in time from one another, and occur instantaneously. Here, the 
collision time is effectively zero, while the average correlation time is the 
mean time between collisions. In the second case, the impact approximation 
holds if the collisions are so weak (although they may last for a significant 
time with their effects even overlapping) that their combined effect on the 
wave is small in the characteristic collision time. 

The quasistatic approximation is based on the assumption that the ex- 
cited atoms are in a relatively constant perturbing field produced by randomly 
situated (slowly moving) perturbers. The resulting electric microfield per- 
turbs the energy levels; this disturbs the emitted line frequency by an amount 
Av, which depends on the strength of the interaction. Thus, the total line 
shape depends on knowledge of the probability distribution for stationary per- 
turbers to provide a given interaction with the atom. The slowly moving ions 
are particularly suited to this treatment; whereas, the faster moving electrons 
may be treated more adequately by the impact approximation, especially in high- 
temperature plasmas. However, as indicated in Section 1.0, INTRODUCTION, it 
has been shown that, for the line wings, the quasistatic approximation has been 
used successfully to describe the broadening because of electrons and ions in 
low-density, low-temperature rf hydrogen and helium plasmas. This was true for 
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the hydrogen case where the empirical correction of Schliiter (ref. 7) was used. 
In the present work, the same type of correction was applied to the quasistatic 
broadening theory of helium. 

Pfennig and Trefftz (ref, k) solved for the helium first-order energy- 
perturbation matrix elements resulting from a static external electric field 
and their associated line strength as a function of the electric field. This 
was done in parabolic coordinates by using essentially the same approximation 
as Foster (ref. l4). They formed the field-dependent off-diagonal matrix ele- 
ments 




+ z. 


n,L - 1 


> 


( 6 ) 


where e is the charge of the electron; z^ and z^ are coordinates of the 
electrons; and n,L,M are quantum numbers. 

These matrix elements form block diagonal submatrices that are asso- 
ciated with each M. The elements on the principle diagonals are the unper- 
turbed energy eigenvalues. For L > 1, the helium eigenstate 4* (n,L,M) may 
be adequately described by an unexcited Is electron which does not contribute 

to H and by an excited electron n£m(& = L,m = M), the eigenfunction of 

l—i " JL ^ 

which is approximated by the hydrogen type. Thus, according to Bethe and 
Salpeter (ref. 15) 



where a is the Bohr radius, 
o 

Pfennig and Trefftz (ref. 9) used the values for the diagonal elements 
H of Martin (ref. l6) and extrapolated them to extend the table. They 
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expressed the Stark-effect eigenfunctions as linear combinations of the unper- 
turbed helium functions i 1 (n,L,M). 

o 


= X) a L ^(E)¥ o (n, L ,M) 


( 8 ) 


Then, they solved the eigenvalue problem 

V H , a , = a £ 

' LL L p Ly p 

with the normalization condition 


(9) 


n-1 

E 


L=M 



1 


( 10 ) 


( U = M,M + l,...n - l) . The index y is related to the parabolic quantum num- 
ber n^, that is, y = n^ + M. Pfennig and Trefftz (ref. 9) tabulated the Stark 
effect, splitting Av (perturbation energy) and relative line intensities . 

The splitting is relative to the unperturbed nD level, and perturbations in 
the 2P levels are neglected. 

The equation for the perturbed relative line intensity of one component 
is given by the equation 


V 4v) 


■ 

^ —CO *— 


4u m„ < e) 


AvJL m ^(E)w(E) dE 


(11) 


where L m ^(E) is the relative line intensity as a function of perturbing field 
strength and W(e) dE is the probability of finding the atom in a perturbing 
field E < E' < E + dE. The integration with the Dirac delta function picks 
out a particular field strength associated with a given wave number. 
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Integration in equation (ll) by the use of the theory of generalized 
functions yields 

I., (Av) = L,. /EL, \W/EL, \ 

My My^ My J \ My J 

where is the electric field that is required to produce a displacement 

Av in the M,y component. A summation of the overall components yields 


dE 


dAv 


My 


( 12 ) 


1<4u) ■ £>( e m,) w ( e m„) 


dE 


dAv 


(13) 


My 


This is the expression used by Pfennig and Trefftz (ref. 4) in their 

calculation. However, this equation is not valid for certain values of Av. 

I dE 


In several of the components , 


dAv 


is singular for certain Av. This is 


My 


because of a reversal or stoppage in the displacement Av with increasing 
field strength. To avoid this mathematical difficulty, equation (ll) was 
folded with the Doppler broadening function. 


D( Av) 


2 2 

M 2 ^ AV 

" 2kT _ _2_ 1 -(Av/o) 

2irkT 6 y/v a 6 


(14) 


One interesting feature of this function is that, in the limit a -+ 0, the 
Doppler function behaves like the delta function. However, in the physical 
situation, the Doppler function has a finite width; thus, the singularity en- 
countered by using a delta function is eliminated. 

The general expression for the folding of two broadening functions is 


P( Av) 



Av')T(Av - Av') d(Av') 


(15) 
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Equations (ll) and (l4) are substituted into equation (15) to yield 


P( Av ) 




(E) 


Av n L,. 

J Mu 


(E)W(E)D(Av - Av') dAv' dE 


(16) 


The integration of equation (l6) with respect to Av' yields 


I ( Av ) 


-r 


w(E) L Mp (E) D 


£av 



dE 


(17) 


where P -*■ I in the notation. Thus, the Doppler function has replaced the 

delta function in the integral. This integral may now be computed by numerical 

means. The desired equation for the triplet lines is obtained by restoring all 

3 

components of the line and by considering the splitting in the 2 P level. 


I( Av) 


V 2J + 1 
2—t Q 
J=0 y 



w(e)l 


My 




(18) 


where Av„ t (E) = Av,, (E) + Av t , Av T = 1.045 for J = 0, Av T = 0.049 for 
MpJ Mp J J J 

J = 1, and Av = 0.029 for J = 2. Because the integrand is defined as zero 

tj 

for E < 0, the lower limit of the integration may be zero. 

If W(E) and D(Av) are normalized and if V'l„ (E) = 1, then 

My 



(19) 


In the tables of Pfennig and Trefftz (ref. 9), ^L Mp (E) = 1 when the 
sharp component of the line (y = 0 term) is omitted. 

The field-strength distribution function W(E) was obtained from the 
H(g,y) tables of Baranger and Mozer (ref. 10 ) for the high-frequency component. 
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In these tables, a two-body correlation approximation was used to account for 
shielding when 3 is the dimensionless field strength E/E . Recently, 

Hooper (ref. 17) has calculated the field-strength distribution function. His 
calculations account for all the higher order correlations; however, for small 
Debye shielding parameters y, the tables of Baranger and Mozer (ref. 10 ) dif- 
fer very slightly from those of Hooper (ref. 17). They tabulated H(g,y) for 
3 from 0.2 to 10 and for y from 0 to 0.8 where y = d/A^. The following . 
equation defines d. 



3/2 3 

d°N = 1 


(20) 


where N is the perturber density and 
length. The normal field strength E q 



is defined by 


is the Debye 


E 

o 



( 21 ) 


Equations (20) and (21) are combined to yield 

E = 3.7^8 x 10 -1 °N 2//3 KV/cm (22) 

o e 

The function W(E) was obtained from the tables of H(3,y) "by the 
relation 


W(E) = jjj- H(3,y) 
o 


( 23 ) 



’’or 


3 > 10 , the asymptotic formula of Baranger and Mozer (ref. 10 ) was used. 


H(S) 


y3 


-7/2 


(o.5 1 +53y + 11.783 1 + ll4.63 "- 5//2 + * * •) 


(24) 


where H « 1.4-963 5/2 + 7.6393 ** + 21.603 11/2 + ••• . 
o 

An empirical correction to the line profile may be made by following the 
reasoning of Griem (ref. 18), Schliiter (ref. 7), and Edmonds, et al. (ref. 8). 
Here, an attempt was made to extend the basic quasistatic Stark broadening 
theory toward the center of the lines by making the effective perturber density 
a function of the wavelength. Effectively, this makes the line core a function 
only of ions , whereas the wings become dependent on both electrons and ions in 
a smooth transition. The empirical correction used is of the form 


N = N 

e 


1 


1.5 + 0.5 



(25) 


and 


z = 



(26) 


where a and b are constants to be chosen, Av < Av is the displacement 

Li 

from the line center for which impact treatment of the electrons is appropriate , 

and Av > Av T is the range of validity of the quasistatic approximation. 

Li 

The value of Av was calculated by the use of the asymptotic behavior 

l_l 

of W(E) and by averaging over Stark components. For the quasistatic approxi- 
mation to be valid, the time between collisions x = r/v must be long compared 
with the perturbation period 1/2tt|Av|c. Here, v is the velocity and r is 
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the distance of the perturbing particle which produces the electric field. 

Thus , 

t >> 1/2tt|Av|c (27) 


Because 



and 


1 2 
2 m e V 



the following condition exists. 


Av >> Av 


L 


1 3kT E_ 

,22 m Av 
4u c e e 


( 28 ) 


This condition is too strict, as a precise calculation of the theory shows. 
It is sufficient to have the following condition (ref. 19). 


Av 


L 


0.2 


1 3kT E 


,22 m Av 
4ir c e e 


(29) 


Values were calculated for E/Av and were averaged over all of the Stark com- 
ponents in the following manner. 


E = 



V E)W (%) 

dE 

dAv 

Mjj 

£V E,W CW) 

dE 

dAv 

Mp 


(30) 


where the factor in the parentheses is a weight factor which weighs the com- 
ponents by the fractional intensity of the component. With the approximation 
s E 

that I — 5a (justified in the linear Stark effect and therefore only 

|dAv My 

-5/2 

slightly dependent on Av) , the asymptotic expression W(E) ~ E 


may be 
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used to obtain the approximate expression for the average field strength, as 
shown in equation (31). 


E 


£V e) V 3/2 


(31) 


Because E is only slightly dependent on the value of Av , E/Av may 
be introduced into the expression for Av . In the numerical calculation, a 
component was omitted if it did not contribute to the line profile in the range 
of interest; for example, the S (p = 0) component was omitted. 

The technique for extending the range, for which the quasistatic theory 
applies, is limited because the profiles are no longer automatically normalized 
if W(E) has been normalized by the use of a constant normal field strength 
E q . To avoid renormalization of the line profile, W(E) was normalized in the 
following manner. Taking the normalization condition 


/; 


H(B)dB = 1 


(32) 


and transforming the variable of integration the following is obtained 



(33) 


where 


3 £ 

9E 


is the absolute value of the Jacobian of the transformation. 
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Therefore , 
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A modified field-strength distribution may be defined as 


W(E) 



E_ 9E 9N 9Av 

E 9N 9z 9Av 9E 
o 


( 35 ) 


(36) 


This differs from the usual definition of W(E) by the Jacobian factor. Here 
the normalization is preserved. 

In calculating the derivatives in the Jacobian factor, equations (22), 
(25), and (26) were used. The dependence of Av on E was found in the tables 
of Pfennig and Trefftz (ref. 9 ) for each Stark component. The derivative was 
calculated numerically by fitting a parabola to three adjacent points about E 
and then by calculating the derivative of the parabola at E. 

In calculating the theoretical line profiles, equation (l 8 ) was pro- 
gramed for computation on a digital computer. The tables of line intensities 
, displacements Av (eigenvalues), and field-strength values E were taken 
from Pfennig and Trefftz (ref. 9)- The tables were interpolated for interme- 
diate values of E in the integration. The field-strength distribution- 
function H($,y) tables of Baranger and Mozer (ref. 10 ) were interpolated 
with the calculated value of y for each $. The resulting distribution was 
then used in the calculation of W(E) according to equation (36), interpolating 


for intermediate values of E. 



3 . 0 EXPERIMENT 


3.1 THE APPARATUS 

3.1.1 The Arc-Jet Facility 

Measurements of the spectral line, continuum intensities, and the line 
profiles were made in the 1.5-megawatt arc-tunnel facility at the Manned 
Spacecraft Center. An arc heater, shown in figure 1, is used to heat various 
gases to high temperatures. In these experiments helium was used. 

The arc heater was designed and built by Electro-Optical Systems, Inc., 
under a NASA contract. The arc heater consists of a tungsten-tipped cathode 
at the upstream end, a restrictor tube with an inside diameter of 2 cm, a di- 
vergent conical nozzle, and a ring of anode pins in the nozzle. All compo- 
nents are water cooled. The restrictor section consists of a number of 
annular-ring, water-cooled segments made of copper. They are insulated from 
one another with boron nitride inserts. Also, the segments are separated by 
0-ring seals. The conical nozzle is segmented in the same manner as the re- 
strictor section. One of the segments in the nozzle consists of a boron 
nitride annular ring and it has 30 radially drilled holes which accommodate 
the copper anode pins. Thus, each of the 30 radial pins forms part of the 
anode. Each pin is connected in series with a 1.4-ohm ballast resistor, which 
insures that the current is divided evenly among the 30 pins. 

Gas is injected approximately tangentially near the cathode through an 
entry block made of boron nitride. There is a secondary-gas inlet approxi- 
mately 12 inches downstream from the cathode. This inlet can be used to in- 
troduce a secondary gas such as oxygen. However, in this experiment, this 
inlet was used to inject helium. A small amount of an inert gas is also 
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Figure 1.- Schematic diagram of arc heater. 
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introduced around the anode pins. This gas cools the pins and the boron ni- 
tride ring which insulates them from the copper segment. 

The power applied to this arc heater is limited to 1500 A and 1 MW. 

The current limitation is because of the cathode. The tungsten tip begins to 
spallate at higher currents. The 1-MW limitation is because of the heat load 
to the restrictor section, which was designed to absorb up to 320 KW of com- 
bined radiative and convective heating. Early in the experiments , it was 
found that 1400 A was the maximum current that could be used reliably because 
of the failure of anode ballast-resistor elements at the points of contact 
with the water-cooling- tube fittings. Later, a modification in the design of 
these fittings eliminated this problem. Because data had been taken at the 
lh00-A condition, it was decided not to increase the current, although running 
at as high a power as possible was desirable because the plasma luminosity in- 
creased with an increase in power. 

The arc heater was designed to run on gas-flow rates of from 2 to 
20 g/sec. A lower flow limit is necessitated by the increase in heating to 
the restrictor section at very low gas-flow rates. It was found that the 
maximum gas enthalpy (calculated by the heat-balance method) could be obtained 
at gas-flow rates of from 2 to 3 g/sec. A total gas-flow rate (not including 
anode-pin bleed) of 2.0 g/sec was chosen. While running on pure helium, a 
flow of 1.5 g/sec was introduced at the cathode, and 0. 5 g/sec was introduced 
in the restrictor section (secondary inlet). An anode-pin bleed of approxi- 
mately 0.5 g/sec of helium was used to externally cool the boron nitride ring. 

Four pressure taps are located on the arc heater to measure pressure. 
One tap is located near the cathode, a second tap is approximately in the 
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center segment of the restrictor, a third tap is in the "throat" or downstream- 
end segment of the restrictor, and the fourth tap is at the nozzle exit. 

The initial arc discharge is begun at a low pressure (approximately 
0.5 torr) and a starting gas-flow rate of approximately 0.5 g/sec. Argon and 
helium have breakdown potentials that are low enough that an open-circuit 
voltage of 1000 volts will initiate the discharge. Once the discharge is es- 
tablished, the main gas flow is turned on and the starting flow is then turned 
off. 

The heated gas is expanded to supersonic velocities through the coni- 
cal nozzle into a vacuum test chamber. The nozzle has a full-divergence 
angle of approximately 56°. The exit diameter is 6 inches. 

The vacuum test chamber and the arc heater are shown in figure 2. 

The chamber is approximately 6 feet in diameter and contains apparatus for 
inserting ablation models and probes into the jet. Mirrors are mounted inside 
to direct light from the jet or ablation models out through windows into light- 
measuring instruments . The chamber is connected to the diffuser section and 
heat exchangers, and then is connected to a four-stage steam-ejector system 
which maintains a vacuum in the system. A large valve, located in the vacuum 
line between the steam ejectors and heat exchanger, is used to control the 
pressure inside the test chamber. The valve is electrically operated from the 
control room so that the chamber can be isolated from the steam ejectors by 
closing the valve, and so that the pressure inside xhe chamber can be regu- 
lated by adjusting the extent to which the valve is opened. In these experi- 
ments , a pressure of k .76 torr was maintained in the test chamber and was 
measured by a radioisotope-type ionization gage (Alphatron). 
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Power to the arc heater was supplied by a 1.5-MW ac and dc power 
supply operated in the dc mode. This power supply will operate in three 
basic dc configurations (table I). Four pairs of modules in parallel can be 
linked to provide the current and voltage combinations shown in table I. 

T'he current is controlled by saturable-core reactors on the ac side of 
the power supply. Rectification is accomplished by the use of three-phase sil- 
icon diode bridge networks . 

The arc heater, test chamber, and heat exchanger are cooled by 500-psi 
deionized water flowing in a closed system. This water is then cooled through 
a water-to-water heat exchanger which is supplied by a secondary open-loop wa- 
ter system. A small amount of water is continually recirculated through 
deionizers. Deionization of the water is necessary to prevent buildup of 
mineral deposits in the arc heater and other water passages and to minimize 
the water conductivity. 

3.1.2 The Facility Instrumentation 
and Measurement Techniques 

The facility instrumentation consists of various types of transducers 
for measuring pressures, temperatures, flow rates, voltages, and so forth. 

The electrical output signals from these instruments are recorded on a 
50-channel analog-to-digital magnetic tape recorder. The magnetic tapes of 
data are then processed on a digital computer in the MSC computation center. 

A sample of the computer output is shown in table II. The time-averaged val- 
ues and minimum and maximum values of each measured quantity are given. The 
variations indicated by the minima and maxima are caused by noise and drift 
(for nominally constant quantities). 
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TABLE II.- SAMPLE DATA FROM DATA REDUCTION COMPUTER PROGRAM 


National Aeronautics and Space Administration 
Manned Spacecraft Center - Houston 
Experimented. Heat Transfer Section 

Statistical Report Number of Samples 1915 
Test Run 423 Run 4230007 Tape 52676 

Time Seconds Date 

Starting 2.6062 12-19-68 

Ending 64.3764 12-19-68 


Measurement Description 

Average 

Units 

Range 

Event 

3.78110 


-89.3786 

94.9964 

Main H^O Temp In 

67.2433 

Deg.F 

53.7282 

82.9703 

Arc Heater H^O Temp Out 

84.2175 

Deg.F 

78.3345 

92.9680 

Delta T 

16.5217 

Deg.F 

16.3581 

16.7089 

Arc Heater H_0 Flow 

C. 

18.3722 

lb/ sec 

18.2266 

18.4924 

Mid- Con Helium Flow 

. 120825-02 

lb/ sec 

. 120611-02 

.121045-02 

Cathode Helium Flow 

.303156-02 

lb/ sec 

. 302850-02 

.303516-02 

Voltage 

475.953 

Volts 

474.277 

478.083 

Current 

1409-71 

Amps 

1403. 55 

1419.88 

Nozzle Pressure 

4.59730 

Torr 

3.93940 

5.35605 

Throat Pressure 

286.135 

Torr 

284.976 

287.074 

Mid-Construction Pressure 

502.873 

Torr 

500.931 

504.665 

Cathode Pressure 

581.800 

Torr 

578.129 

585.840 

Chamber Alphatron 

4.82645 

Torr 

4.78047 

4.89595 

Nozzle Alphatron 

3.95931 

Torr 

3 . 9 l 8 l 6 

4.00741 

Intensity 

-7 . 4l42l 

MV 

„44.3476 

.607147-01 

Abel Scan Position 

5.344o8 

MV 

-34.2110 

36.7661 

Abel Scan Command 

57-7177 

MV 

-.833321-02 

94.2833 
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The measurements of the spectral lines were made at one set of run 
conditions; therefore, the values shown are considered standard and typical. 

The instruments used to measure facility parameters are described as 
follows. Copper-constantan thermocotiples are used to measure water tempera- 
tures. The temperature increase in the water which flows through the arc 
heater is measured by a thermopile, the output of which is a millivolt signal 
that is approximately proportional to the temperature difference between the 
water inlet and the water outlet. 

Water- flow rates are measured by turbine-type flowmeters, the output 
of which is a frequency that is proportional to the flow rate. This sinus- 
oided signal is converted to a dc signal by a frequency-to-analog converter. 
Gas-flow rates of the primary gases are measured by the sonic-orifice tech- 
nique. The pressure upstream of a calibrated orifice is measured with a 
strain-gage-type pressure transducer. The flow rate is then calculated by the 
use of the sonic-orifice equation 



( 37 ) 


where m is the mass-flow rate, A is the effective orifice area, P is 

eix 

the upstream pressure, R is the perfect gas constant, T is the temperature, 
and y is "the ratio of specific heats of the gas. For this formula to be valid, 
the upstream pressure must exceed the downstream pressure at least by about a 
factor of two. 

The current through the arc heater was determined by measuring the 
potential drop across a calibrated resistor which was in series with the arc 
heater. This voltage was transduced by an isolating device, the dc output of 


27 


which is proportional to the dc input. The potential drop across the arc 
heater was measured by a similar transducer. Isolation of the recorder from 
the arc heater power supply was necessary for safety reasons and to protect the 
recorder from possible large-current ground loops. 

The nozzle-exit pressure and the test-chamber pressure were measured by 
an Alphatron ionization gage. The other pressures in the system were measured 
by strain-gage-type pressure transducers. 

A measure of the net energy introduced into the gas by the arc heater 
is the heat-balance enthalpy. It is defined by the equation 

IV - C m AT 

h = — ( 38 ) 
m 

g 

where I is the current across the arc heater, V is the voltage across the 
arc heater, is the specific heat of water, AT is the temperature in- 

crease of the cooling water through the arc heater, m is gas-flow rate, and 

g 

m w is the cooling-water flow rate through the arc heater. Strictly speaking, 
this is not the actual thermodynamic enthalpy, because it is a quantity aver- 
aged over the length of the arc heater and the pressure is not constant over 

the length of the arc heater. In these experiments h = 1.8l 4 0.05 x 
12 

10 erg/g. The standard error is given here along with the calculated values. 

Another measure of the stagnation enthalpy is given by the sonic- 
throat equation (ref. 20 ). 



m 


( 39 ) 
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This equation is based on the assumption that the gas behaves as a perfect 

2 

gas. P q is standard atmospheric pressure, 2115 lb/ft ; A* is the area of 

2 

the sonic-nozzle throat, ft ; and P is the stagnation pressure, atmospheres. 

The two numbers in parentheses are conversion factors. In these experiments 

12 

h s =1.95 * 0.09 x 10 erg/g. This technique for determining the stagnation 
enthalpy is only an approximation because of the following factors. The gas 
is not strictly perfect, because there is some ionization, Wo true stagnation 
pressure exists , because there is a pressure gradient along the length of the 
restrictor. The pressure measured at the cathode was used in the sonic-throat 
equation; it is considered that this pressure is most representative of a 
stagnation pressure, because the gas velocity is lowest at that location. 

3.1.3 The Spectroscopic Instrumentation 
The spectroscopic data were obtained by the use of a Jarrel-Ash 3.4-m 
Ebert stigmatic spectrograph. This instrument has both photographic and photo- 
electric capability. A grating having 590 lines/mm and an effective aperture 
of llU mm was used, yielding a total of 67 260 rulings. The spectrograph has 

O 

a reciprocal dispersion of 5-05 A/mm in the first order with this grating. 
Photographic plates were used to identify the spectra and for qualitative 
evaluation of the spectra. 

A sine bar-grating drive rotates the grating so that the spectral line 
scans across the exit slit linearly in time or with shaft rotation when a 
constant-speed motor drives the sine-bar lead screw. Various combinations of 
entrance-slit and exit-slit widths were used during these experiments. The 
choice depended on the type of measurement and on the width of the spectral 
lines being measured. It was desirable to use as wide a slit as possible, 
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narrow enough to sufficiently resolve the line in question. It was fortuitous 
that the weaker lines were broader, so that wider slits could be used to in- 
crease signal strength. 

Photoelectric measurements were made using an EMI-6255B photomulti- 
plier tube that has an S-13 response. The anode was connected through a 1M9 
resistor to ground. 

The resultant signal was connected to the input of the phase-sensitive 
lock-in amplifier. The light going to the spectrograph was interrupted at a 
rate of 10 800 Hz by a light chopper that was mounted next to the entrance 
slit. The light chopper had an output signal that was used as a reference by 
the amplifier. Positioning the chopper in front of the entrance slit was very 
critical, because openings in the blades were small. There were 180 openings 
near the periphery of the 9-inch-diameter chopper blade. Vibration induced 
into the spectrograph when the chopper was mounted directly to the optical 
bench was serious enough to cause fluctuations in the output signal when using 
narrow slits. To eliminate this problem, the chopper was mounted on a support 
which was isolated from the spectrograph. A schematic of the arc tunnel, 
spectrometer, and optical path is shown in figure 3- 

The optical path passed from the arc jet F c through two mirrors, 

and M^. These two mirrors were needed so that a region 5-inches from the noz- 
zle exit plane could be viewed perpendicularly to the axis of the jet. The 
window in the tunnel had a diameter of 7 inches , and its center was 9 inches 
from the nozzle exit plane. The window was made of fused silica. 

After the light passed from the tunnel, the image was rotated 90° in 
a plane perpendicular to the optical axis by a pair of mirrors . At the same 
time, the optical axis was directed to the scanning mirror M_. Rotation of 



To heat 
exchanger 
and steam 
ejectors 



Figure 3.- Arc tunnel, showing the optical arrangement. 
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the image was necessary so that the image of the horizontal jet could be made 
to lie along the length of a vertical entrance slit. 

The scanning mirror was driven by a stepping motor, and its angular 
rotation was sensed by a single-turn resistive-film potentiometer which had a 
linearity of ±0.01 percent. The angle of rotation of the mirror required to 
scan from one side of the jet to the other side (a distance of lU cm) was 2.6°. 
Therefore, the accuracy in mirror position was reduced to ±1.^ percent of full 
deflection, neglecting any other sources of error. 

From the scan mirror, the light was directed through a 70-mm-diameter 
collecting lens, L^, then to a 30-min-diameter focusing lens, L^, then to the 
chopper and entrance slit of the spectrometer. All lenses were made of fused 
silica to provide transmission in the near ultraviolet. All mirrors were 
first-surface mirrors deposited with pure aluminum and were without overcoat. 

It was found that the mirrors inside the test chamber had to be replaced oc- 
casionally, because a film deposited on their surfaces. This film was prob- 
ably a product of ablation which had condensed onto surfaces of the test 
chamber during previous runs during which ablation materials were tested. 
However, it should be pointed out that no evidence of ablation products ap- 
peared in the long-exposure spectrograms. The only identifiable impurities 

that were detected in the spectra were a weak OH band in the spectral range of 
o + o 

3063 to 3^-00 A and a very weak bandhead near 3900 A. This probably came 
from water vapor which had diffused into the test chamber from the steam ejec- 
tors and small air leaks. 

Reference sources and calibration sources were placed at the focal 
plane F . Likewise, an iron arc was located there when the optics were being 

n 

alined. A ribbon- filament lamp was also placed at the focal plane F c , in the 
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chamber for use as a calibration source to determine the spectral transmission 
of the optical path. To determine the largest slit height usable while cali- 
brating with the source at F , the photomultiplier output was measured com- 
pared with the slit height. The results were plotted, and the curve was 
linear up to a height of approximately 7 mm. Therefore, a 7-mm slit height 
was used when calibrating the optical path. 

To verify that the transmission function of the optics would remain 
the same for any portion of the jet in a plane perpendicular to the axis of 
the jet, the lamp was centered approximately 10 cm above, 10 cm below and on 
the jet axis. The variation in the signal from one position to another was 
less than ±1 percent. The radius of the luminous jet was 2.75 inches (6.99 cm) 
at the position 5 inches from the nozzle exit plane. This was determined by 
inserting a water-cooled probe into the jet at the lower visible boundary, 
then at the upper boundary, and taking the difference in probe position. The 
probe was positioned electromechanically by a servo system. 

To obtain Abel-inverted line profiles, the scan mirror and the grating 
were positioned sequentially by a programable stepping-motor indexer. This 
device drove both the grating-scan stepping motor and the position-scan step- 
ping motor. The motor shaft rotated 1.8° per step. The motors were connected 
to the grating drive and the scan mirror through reduction gears to obtain the 
very small rotations necessary for scanning. A magnetic clutch was used in 
the grating-drive linkage so that the grating could be rotated manually from 
one line to another without having to disconnect the motor. A drawing of the 
(Abel) position-scan-drive assembly is shown in figure 4. 

The position- (or Abel) scan motor was stepped a preset number, N 
steps, M times. After each N steps, the motor stopped for a preset time, t. 







The duration of time was determined by the time constant of the amplifier. 

A longer time constant was used when scanning weak lines. M of the stops com- 
pletes one scan of the jet. Then, the grating was incremented K steps, cor- 
responding to the desired wavelength increment. The wavelength increment was 
determined by the breadth of the line and the resolution desired. This whole 
sequence was repeated L number of times. L was determined by the number of 
increments required to scan across the line, or part of a line in the case 
that a single line was scanned piecewise. A piecewise-scan technique was used 
to obtain a different sensitivity range of the amplifier for different por- 
tions of the line. While the motors are stopped at each position during a 
scan of the jet, the indexer outputs a dc signal to one channel of the re- 
corder. This signaled the computer when to accept the intensity and position 
data which were recorded in other channels. The position was sensed by the 
potentiometer previously mentioned. Each fraction of the potentiometer re- 
sistance formed one side of a Wheatstone bridge, such that the output of the 
network was proportional to the position (fig. 5). 

The scan-mirror assembly was calibrated by focusing an arc, located at 

the focal plane, F^, onto focal plane F c> The distance the image traveled 

in a given number of steps was measured several times. The average resulted 

in a distance of 0.0398 ± 0.000^ cm/step. The grating drive was calibrated by 

rotating the shaft a certain number of turns (or steps), and recording the 

corresponding wavelength difference. This resulted in a calibration of 

o 

0.0025 1 1 percent A/step. 
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3.2 THE SPECTROSCOPIC MEASUREMENTS 

The region of the jet in which the spectral measurements were made is 
shown in figure 6; the rectangular overlay denotes the region scanned by the 
entrance slit for the Abel-inversion measurements. The width of the rectangle 
is the height of the slit image projected onto the jet. This region is down- 
stream of an expansion shock wave, which is also shown on the overlay. The 
outer edge of the jet cannot be seen in this photograph because the exposure 
was reduced to reveal the shock wave more clearly. 

3.2.1 Line Intensities 

The total-intensity jet profiles of several lines were measured by the 
use of the Abel-scan mechanism. These measurements were made for a wide exit 
slit (205y ) and a somewhat narrower (50y) entrance slit. This facilitated 
measuring the total intensity at the peak wavelength (calibrations are given 
in appendix A). 

A scan across the line was made to determine the extent of the line in 
wavelength and to determine the constancy of the background continuum in the 
neighborhood of the line. Then, the line was centered on the exit slit and an 
Abel-scan across the jet was made. Another scan was made of the background 
continuum a few angstroms from the line. After performing the numerical Abel- 
inversions of these data, the Abel-inverted-background continuum signal was 
subtracted from the Abel-inverted total line-intensity signal. The net signal 
was converted to the total line intensity in proper units by the use of equa- 
tion (Al6) . 
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Figure 6.- Photograph of portion of jet which is scanned by the spectro 


The conical expansion shock wave is denoted, and it makes an angle of 


approximately 30° with the jet axis. The outer edge of the jet is 


not visible here because the exposure was reduced to enhance the 


appearance of the shock wave in the original photograph 
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3.2.2 Line Profiles 

The shapes of a number of line profiles were measured with the center 
line of the arc jet focused on the entrance slit. This measurement yielded an 
average profile over a diameter of the jet. Two techniques were used to scan 
the line. The first method involved a scan of the entire line at a constant 
stepping rate and on a single amplifier-sensitivity setting. The second meth- 
od involved a scan of the line piecewise, changing the amplifier sensitivity 
to take full advantage of the amplifier accuracy on each range. The stepping 
rate and amplifier time constant were adjusted to obtain the best signal-to- 
noise ratio of the output consistent with a reasonable rate of acquiring the 
data. The scan was continued until a minimum or a constant signal was ob- 
tained. This minimum constant signal was taken to be the background- continuum 
signal level. The continuum signal was subtracted from the scan signal to ob- 
tain the true spectral profile of the lines. In the case of weak lines, this 
was difficult to do accurately because the signal becomes less than the noise. 
Therefore, in the case of weak lines, the accuracy in the far wings suffers. 

In addition, overlapping of lines caused difficulties in some cases. 

3.2.3 Abel-Inverted Line Profiles 

To obtain profiles that represented the line shape in a homogeneous 
plasma (that is, a profile which is not an average over inhomogeneous layers), 
it was necessary to Abel-invert the profiles. The apparatus and the sequence 
of the scans are described in Section 3.1.3, The Spectroscopic Instrumentation. 

A sample of the data is shown in figure 7; the intensity, position, 
and command signal are shown. The units are in millivolts compared with time 
(sec). These raw data were reduced by computer according to the following 
scheme. The data just mentioned were read off the magnetic tape in digital 



Position, Intensity, 

Command signal , mV arbitrary units arbitrary units 
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150 160 170 180 190 200 210 220 230 


Elapsed time, sec 


Figure 7 • - Sample Abel profile data showing scans of the jet at 3 wavelengths 
in a line profile out of a total of approximately 100 required to make a 
complete Abel-inverted line-profile measurement. Fifty-eight steps are 
in each Abel scan. The command signal >90 mV denotes a point in the scan. 



form. Only the data that were flagged by the command channel were sampled. 

In addition, the position was required to be within the jet radius R. A 
counter kept track of the number of samples in a scan of the jet so that this 
process was repeated for each scan. By the use of calibration factors, the 
raw position data were converted to centimeters and the intensities and cor- 
responding positions were stored on tape for each scan. 

To reduce the uncertainty in determining the center of the jet, the 
center position was found by calculating the centroid 


y = XXXi) (4o) 

° XX 

where l(y ) is the measured intensity at the corresponding position y^. 

The corrected positions were then yj = y_^ - y . The light scattered and re- 
flected by the walls of the chamber was accounted for by scanning beyond the 
jet radius, R, approximately 1 cm above and below the jet, and averaging 
those measured intensities which were outside the jet diameter. This average 
was subtracted from the intensity data. 

The data were then Abel-inverted according to the method described in 
appendix B. 

3.2.4 Continuum Measurements 

o 

The series limit continuum was scanned from 3500 to 3200 A. This was 
done by the use of an exit slit width of T5y- The absolute calibration is 
described in appendix A. One difficulty in making this measurement was the 



4i 

o 

rather weak signal in the range of 3300 to 3100 A and an impurity band OH 
(A E - X it), which has a bandhead at 3063 A and is degraded toward the red 
(fig. 8), 

The Abel-inverted absolute-intensity profiles of the jet were measured 
in the series-limit continuum and in the continuum near each of the Abel- 
scanned lines, as was mentioned in Section 3.2.1, Line Intensities. These 
latter measurements were made with an entrance slit of 50 p and an exit slit 
width of 205 m- Inasmuch as the intensity of the continuum was weak, it was 
necessary to scan slowly. This required averaging times of 1 to 3 seconds per 
position to smooth the data. 


3.2.5 Instrument Profiles 

To account for the instrument broadening in the measured profiles , it 
was necessary to measure the instrument profile for each configuration that 
was used in the experiments. It was found that the instrument profile need 
not be taken into consideration in a number of cases, because many of the 
lines were much broader than was the instrument profile. The instrument pro- 
file can be measured by profiling a spectral line whose width is significantly 
less than that of the instrument profile. The profile which is measured is 
the folding of the line profile with the instrument profile 

1(A) = f P(X')S(X' - x) dX ' (la) 

Jo 

where P(l) is the instrument profile and S(A) is the profile of the source 
line. Because S(X) is much narrower than P(X), it acts like the Dirac 


delta function in the integration, so that when carrying out the integration 
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Figure 8.- Photographic spectra of arc jet. Parts (a) and (d) are spectra 

of a mixture of 97 percent He and 3 percent H^. Parts (b), (c), (e), and 

(f) are pure helium spectra, made with differing exposures and slit widths 

to show both faint and strong lines, respectively. Parts (d), (e), and 

(f), are the longer wavelength extensions of parts (a), (b), and (c), 

respectively. The sharp, faint lines adjacent to the primary spectra 

are made with a He low-pressure discharge for reference. The "Forbidden" 

3 

designation refers to the n P Stark component of the corresponding 

3 3 

2 P - n D lines. 
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I(x) = P( A) is obtained. This justifies the assumption that this measured 

profile can be used as the instrument profile. 

The line used in these experiments for measuring the instrument pro- 

° pgg 

file was the 4358-A line of a low pressure Hg radiofrequency discharge. 

The line from this source is one of the narrowest available, having a width of 

o 

only approximately 0.004 A. The lamp was excited by a 60-watt transmitter 
operating continuously at a frequency of approximately 28 MHz. 

To verify that the spectrometer was correctly alined and that the 
source was indeed narrow, the widths of the measured profiles were compared 
with those calculated theoretically. The theoretical curves of Rolf Brahde 
(ref. 21 ) were used. These curves are the result of folding the entrance and 
exit slits with the diffraction profile of a rectangular opening. The grating 
was the limiting optic in this case. 

The measured instrument full half widths are compared with the theo- 
retical in table III. The measured values are consistent with the theoretical, 
except in the case of the 25 m slit. Good alignment was essential for obtaining 
the narrowest and most symmetrical profiles, and the focus was also critical. 
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TABLE III.- COMPARISON OF THE EXPERIMENTAL AND THEORETICAL 
INSTRUMENT PROFILE FULL WIDTHS AT HALF HEIGHT 


Entrance slit , y 

Exit slit, y 

O 

Measured, A 

O 

Calculated, A 

(a) 

Cb) 

10 

10 

— 

0.051 

0.075 

25 

10 

0.15 

0.126 

0. 14 

25 

25 

0.25 

0.126 

0.15 

75 

75 

0.44 

0.379 

0.43 



A = 4358 A 




Calculated by the use of simple rectangular-slit approximation not taking 
diffraction into account. 

Calculated taking diffraction into account (ref. 21). 
















4.0 DISCUSSION OF RESULTS AND COMPARISON 


OF EXPERIMENT WITH THEORY 

4.1 LINE INTENSITIES AND TEMPERATURE PROFILES 

4.1.1 Determination of Temperature Profile 
from Line Intensities 

To determine the electron temperature distribution in the jet as a 
function of radius, eight narrow helium lines were scanned in wavelength us- 
ing a wide exit slit. These measurements were made at a single jet position 
through the jet center line. From these measurements, an effective spectral 
slit width w was determined by use of the relation (appendix A) 



(42) 


Use of a wide exit slit and this calibration facilitated the use of a single 
peak-wavelength measurement at each position. The character of the back- 
ground continuum was also determined, and it was found that the continuum 
essentially was constant in the neighborhood of each line. 

The jet was scanned perpendicular to the axis to obtain the radial- 
intensity profile of each line by Abel-inversion. Similarly, the radial 
profile of the background continuum of each line was obtained. The resulting 
intensities were calibrated by use of the methods described in Section 3.2.1, 
Line Intensities and appendix A. The temperatures at 10 radial positions 
were determined by making a least-squares fit of the data for each radius to 
the logarithm of equation (l). The slopes of these straight lines equal l/KT. 
The radial temperature distribution is shown in figure 9- The error bars 
denote the standard deviations, which at the center line, amount to an error 



Electron temperature 



Normalized radius, r/R 

Figure 9.- Electron temperature, determined from Boltzmann plots corresponding 
to each radius r/R. Error bars denote the standard deviation. 



of 17.3 percent. The determination of the standard deviation was made using 
standard error-analysis techniques such as those found in Barford (ref. 22). 
These errors are "based on the scatter in the final results of the calibrated 
intensities, and, therefore, they include all sources of random errors. It 
is assumed that any systematic errors are small in comparison. Within the 
limits of experimental error, there was no curvature to the plot of log 
(IX /gf) compared with E , where I is the relative total line intensity, • 
g is the statistical weight of the lower energy level, f is the absorption- 
oscillator strength, and E^ is the upper energy level of the transition. 

The fact that there was no curvature is consistent with the assumption of 
local thermodynamic equilibrium of the electrons in the energy levels involved. 
The lines involved in the temperature determination and the associated 
absorption-oscillator strengths f, statistical weights of the lower energy- 
level g , and the energy in cm ^ of the upper energy level are shown in 

Table IV. The oscillator strengths given here are from Green, et al. (ref. 23), 
and they are the most recent and most extensive tables available. The oscilla- 
tor strengths of Trefftz, et al. (ref. 24) are consistent with these, but 
their tables are less extensive. 

4.1.2 Optical Thinness Verification 

To verify that the optical thinness assumption is valid, an approx- 

3 3 

imate value of the optical depth was calculated for the 2 P - 6 D line and 
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3 3 

for the 2 S - UP line, Griem (ref. 13) gives the following expression 
for the optical depth for a cylindrical source. 


T 




1/2 


k 1 dx = 




rdr 



(43) 


The expression is the integral of the absorption coefficient along the path 
of observation. e(r) is the emission coefficient, determined from an Abel- 
inversion of the intensity scan of the line, and I is the Planck function 

V r> = [ exp Mitt) - x ] (W 


The temperature distribution T(r), which is shown in figure 9 , was fit by 
the least-square technique to a power series. 

7 

t(*o = E V n (^5) 

n=0 


This approximate temperature-distribution function and the function obtained 
from the Abel-inversion calculation were substituted into equation (43). A 
numerical evaluation of the integral was performed. The resulting optical 
depths were 3 x 10 ^ for the 2^P = 6^D line and 2.5 x 10 - ^ for the 2^S - 4^P 
line. These values are much smaller than unity. Therefore, the assumption 
that the jet is optically thin is justified for these lines. Because these 
are low quantum number lines, the jet was assumed to be optically thin for 
all lines that were measured. 
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TABLE IV.- LINES AND CORRESPONDING ATOMIC PARAMETERS 


o 

Wavelength, A 

" 

Transition 

g 

f 

(a) 

BBI 

3187-7 

2 3 S - 4 3 p 

3 

0.0271 

191217 

3613.6 

2 1 S - 5 1 P 

1 

.0231 

193942 

3652. 

2 3 P - 8 3 S 

9 

.00063 

196461 

3732.8 

2 3 P - 7 3 S 

9 

.00102 

195868 

3964.7 

2 1 S - 4 1 ? 

1 

.0506 

191493 

4120.9 

3 3 

2 J P - 5 S 

9 

.00373 

193347 

4437.6 

2 1 P - 5 X S 

3 

.00317 

193664 

4713. 

2 3 P - 4 3 s 



9 

.0104 

190298 


a Oscillator strengths are from velocity calculations by Green, 
Johnson, and Kolchin (ref. 23). 

°Energy levels are from Martin (ref. 1 6 ). 
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4.2 ABSOLUTE CONTINUUM MEASUREMENTS AND THE 
ELECTRON-DENSITY PROFILE 

The electron density profile of the jet was found by solving equa- 
2 

tion (2) for N g . The assumption was made that N g = N^. The following 

data were used in calculating the continuum emission coefficients. The 

photoionization cross sections were those calculated by Anderson and Griem 

(ref. ll), and the bound- free averaged Gaunt factors g^(A) are those 

calculated by Karzas and Latter (ref. 12). The free-free average Gaunt 

factor was set equal to a constant g (A,T) = 1.1. The results of the 

calculations are plotted in figure 10 for T = J800° K and 
15 -3 

N = 1.38 x io cm . Measured values of the emission coefficient were 
e 

combined with equation (2), giving the relationship 

( 46 ) 

2 

where e (A,T) is the measured emission coefficient and [e(A,N ,T)/N ] 

m ’ e’ e c 

is the theoretical value obtained from equation (2). Equation (2) was 

2 

programed for computation on a digital computer. Although e(A,N ,T)/N 

e e 

is a function of N g , its dependence on N g is so small at most of the 
wavelengths of interest that it is negligible over the range of densities 
involved. Therefore, only an order of magnitude guess was needed for the 
first iteration. 

These emission coefficients were calculated as a function of radius 

at the measured wavelengths by the use of the temperature distribution shown 

l4 _3 

in figure 9, and an initial electron density of 6 x 10 cm . Then, from 
equation ( 46 ), electron density distributions were calculated for each 




? -20 

£ (\, t)/|\T x 10 erg/sec 



Wavelength , m- 

2 

Figure 10.” Calculated continuum-emission coefficients divided by N . G 

e 

notation + beside the term designations implies that the curve includes 
the continuum resulting from transitions into that level plus all other 


energy levels above it . 



wavelength and were normalized to unity at r - 0. The measurements for each 
wavelength were then averaged for each value of r to obtain a first iter- 
ation radial-electron density distribution. 

Using this distribution for N g and T(r) distribution, a second 
£ (A,N e ,T) 

set of — p— was calculated. From these, the average radial-electron 

N 

e 

density distribution was determined as in the first iteration. 

The normalized profile was calculated from Abel-scans of the continuum 

. o 

near the wavelengths given in table IV, plus the wavelength 3400 A, which is 
3 3 

in the 2 P - n D series limit continuum. The center line or peak value was 
determined from Abel- inverted measurements at 3184, 3400, 4l2T , 4442, and 

o o 

4718 A plus noninverted measurements at 3250, 3280, and 3475 A. Here, the 

center-line emission coefficient was obtained from the center-line intensity 

by dividing by the mean effective thickness t of the plasma, as shown in 

the expression e(o) ~ l(0)/t. The value t = 3.80 cm was obtained from the 

average value of l(0)/e(0) for each Abel-inverted continuum measurement. 

15 -3 

The center-line electron density was found to be 1.38 x 10 cm ±12.5 per- 
cent. The electron-density profile is shown in figure 11. Error bars are 
plotted at r/R equal to 0, 0.3, 0.6, and 0.9* 

A portion of the series limit continuum and the merging lines are 
plotted in figure 12. The intensity measurement is for y = 0 (along a 
diameter), and the scale is in arbitrary units, with no account made for the 
change in sensitivity of the spectrometer with change in wavelength. 

4.3 REMARKS ON THE ABEL- INVERSION TECHNIQUE 

The technique for scanning the jet and spectrum seemed well suited 


for these experiments because the variable scan rates permitted scanning 



Electron density , N x 10 cm 
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different parts of the line profile at rates which give the optimum signal 
averaging and resolution. For example, in the core of the line where the 
signal is strongest hut higher resolution is required, one could scan rapidly 
with small AX. However, out on the wings where there is less need for high 
resolution, one could scan more slowly a fewer number of points with larger 
AX. This technique permitted a fairly precise knowledge or measure of the 
positions in the scan for each wavelength increment. The estimated error in 
position was about 2 percent of the radius R. The curve fit of the data in 
the numerical Abel-inversion yielded a standard deviation of less than 

1 percent of the center-line intensity near the line core, to approximately 

2 percent on line wings and in the continuum. The fractional error increased 
near the edge of the jet where the intensity is small. The major sources of 
error are related to the jet itself, rather than to the measurement or 
numerical techniques. First, the jet is not strictly symmetrical and de- 
viations from cylindrical symmetry are largest near the edge of the jet. 

There was sometimes a factor of two difference in the intensity near one 
edge compared with the other edge. The technique of fitting the data to 
symmetrical curves effectively averaged the data. In the case of the Abel- 
inverted line profiles, the antisymmetrical component was calculated, but 
there were large deviations in the curve fit caused by a poor choice of 
functions for the curve fit. For that reason, these data are not presented. 
The other major source of error in the radial-intensity profiles stems from 
uncertainty in the overall jet diameter. It was difficult to determine 
accurately where the intensity of the jet decreased to zero because of 
background-scattered radiation. Both visually measuring the luminous 
diameter of the jet and a stagnation-pressure probe survey (fig. 13) yielded 




'ajnssojd uo!ibu6e:|$ 


probe. The center of the distribution is approximately 0.25 inch below the 
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a diameter of l4.0 ± 1.2 cm; hence, the standard error in the radius R was 
± 4 percent. This uncertainty in the radius yielded a similar uncertainty 
in the Abel-inverted intensity W(r). The combination of the errors just 
mentioned yields a net uncertainty of about 5 percent near r = 0, and this 
uncertainty increases to about a factor of two for weak parts of the spectrum 
near the edge of the jet. 

4.4 HELIUM LINE-BROADENING AND 
QUASISTATIC-THEORY COMPARISON 

Experimental line profiles were measured by the use of two techniques. 
In the first technique, the profiles were averaged over a diameter of the jet 
by a line-of-sight measurement. In the second method, the profiles were 
Abel- inverted by scanning both in wavelength and position. 

. 4.4.1 Averaged (Non- Abel-Inverted) Line Profiles 
The averaged profiles are not suitable for precise comparison with 
the theory because they are averages and an accurate comparison would re- 
quire that the theoretical profiles also be averaged over the temperature and 
electron-density distributions along a diameter of the jet. Such theoretical 
calculations would be very lengthy even on a computer, and the results prob- 
ably would not be very helpful in determining the validity of the theory in 
detail. However, for a rough comparison, the experimental average profiles 
and the theoretically calculated profiles are presented in figures l4 to 18 
for the triplet diffuse lines with principle quantum numbers n = 4 to n = 8. 
In these comparisons, the mean temperature 8750° K and a perturber density 
1.58 x 10 cm (which is two times the root-mean-square electron density) 
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Figure 15 . - Concluded 
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Figure IT-- Concluded. 
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were used. The root-mean- square value was chosen because the line intensity 
depends on the square of the electron density because N g = N. . 

In the theoretical calculations, no variation in the perturber density 
with wave number was used, because such a correction is much smaller than the 
accuracy of the overall comparison. The higher lines have uncertain back- 
ground, as was mentioned previously. In addition, the higher lines overlap, 
making the comparisons in the line wings less accurate. For the case n = 4, 
the theoretical profile was folded with the instrument profile. This was not 
necessary for the higher lines, because their half widths exceeded the in- 
strument half width by at least a factor of ten. The profiles for the above 
lines , plus those for n = 9 and n = 10 are shown plotted together in 
figure 19. The profiles in figures l4 to 18 are normalized to the total 
intensity, as shown in the following equation. 

OO 

/ l(Av)dAv = 1 (47 ) 

_co 

Because the profiles do not extend to plus and minus infinity, a correction 

term was added to the integral. It was assumed that the line wings have 
- 5 /2 

a Av dependence; therefore, that part of the wing omitted in an inte- 

gration over a finite interval approximately equals 2/3Av f ,l(Av f ) , where Av f 
is the finite integration limit. The profiles in figure 19 are normalized 
to the peak intensity. 

Allowance for the background continuum and overlapping was made by 
subtracting an estimated background from the profiles prior to normalization. 
However, this was a rather crude approximation for the higher lines. For 
nonoverlapping lines, it would be possible to fit the far wings to a 




(AV) | 
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-5/2 power-law curve and thus determine the amount of background, if one assumes 
the -5/2 power holds for the far wings. If one does not want to make this 
assumption, then the only course is to try to actually measure the background. 

The theoretical profiles for the n = 4 and n = 5 lines exhibit 

strong peaks in their violet wings fairly near Av = 0. In the case of 

n = 4 (fig. l4), the peak at Av = 7 cm 1 is the F q 1 2 ®’ tark - components 

(y = 3). Unlike the experimental profile, the group of F q p p components 

in the theory is well isolated. This indicates that another broadening 

mechanism is needed to establish agreement. In the case of n = 5 (fig- 15), 

the violet wing exhibits two isolated peaks: one at Av = 5 cm ^ and another 

at Av = 11 cm 1 , corresponding to the F Q ^ ^(y = 3) and G Q ^ 2 (y = 4) 

Stark components, respectively. Again, there is a large deviation from the 

experimental profile, indicating the inadequacy of the quasistatic theory 

here. There, is poor agreement in the violet wing because of overlapping 
o 

with the 4009 A line. There is fair agreement in most of the red wing for 
both n = 4 and n = 5- The S(y = 0) and P(y = l) components are not 
present, because they are off the scale of the figures. The agreement of 
the quasistatic profile with the experimental profile for n = 6 (fig. l6) 
is rather poor for two reasons. The violet wing of the measured profile 
has some background interference from an unidentified source yielding a peak 
at approximately Av = 50 cm \ The red wing agreement is disturbed by the 
presence of the forbidden P component. 

For n = 7 and n = 8 (figs. 17 and 18, respectively), the agreement 
in both the red and violet far wings is reasonably good, considering that 
these profiles are averages. The profiles disagree near forbidden components 


and in the core. 
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4.4.2 Abel-Inverted Line Profiles 
Abel-inverted line profiles were measured for three of the triplet 
diffuse lines of helium having n = 6, n = 7, and n = 8. These lines were 
chosen as a compromise. The higher lines had significant overlapping and 
were weak. The lower lines were so narrow as to require narrow exit slits 
to avoid significant instrument broadening, reducing the signal strength. In 
addition, the quasistatic theory does not apply to the low lines until ex- . 
tremely far out on the wings. The lines were scanned as described in Sec- 
tion 3.2.3, Abel-Inverted Line Profiles, and Abel-inversions were calculated 
at each Av, beginning with Av = 0 and progressing toward higher | Av | . 

The scan positions and intensities were then symmetrized (according 
to the scheme described in appendix B) by the use of a linear interpolation 
to obtain intensities at positive positions equal to the absolute value of 
the negative positions. The symmetrization reduced the error in the Abel- 
inversion caused by uncertainty in the center of the jet (ref. 25). Once the 
data were symmetrized, a least squares curve fit was performed on the sym- 
metrical component. It was found that nine terms in the expansion gave the 
best overall fit. The functions had the form (R - y ) , where k took 
all integer values from zero to eight. A typical example of the measured 
data and the fitted curve are shown in figure 20. The corresponding radial- 
intensity distribution, numerically obtained from the curve-fit expansion 

coefficients and the analytic form of the inverted functions (which is 
2 2 k— 1 /2 

A^(R - r") _ ' , where A is a constant resulting from Abel-inversion of 

2 2 k 

(R 1 " - y ) are shown in figure 21. Note that this function is singular 
at r = R for k - 0; hence, the computed radial-intensity distribution is 
invalid in the neighborhood of r = R. 
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After the radial intensity distributions N(r) were calculated for 
each AX , the data were merged to obtain the line profiles at each radius 
r = 0, 0.1R, 0.2R, . . . , 0.8R. 

The profiles including the background continuum and so forth are 
plotted in figures 22 to 2^. Note the wavelength shift in the peaks of the 
profiles with radius. 

The background intensity was determined from other measurements and. 
was assumed to be constant in the neighborhood of the lines. The center-line 
values were subtracted from the r = 0 Abel-inverted line profiles prior 
to plotting in figures 25 to 26. Along with these profiles, several quasi- 
static theoretical profiles are plotted. In these curves, the scale has 
been adjusted so that the maximum intensity is at Av = 0 . 

The profiles for n = 6 are plotted in figure 25. The agreement is 

3 3 

only fair here because of the presence of the 2 P - 6 P component in the red 

wing, and because of the uncertainty in background. The discrepancy in the 

line core where the quasistatic theory does not apply is quite apparent. 

3 

However, even here the influence of the splitting of the 2 D state by the 
Stark effect can be seen. Figure 26 is an expanded plot of the n = 6 line 
core and is presented merely to show the gross correspondence of the peaks in 
the quasistatic theory with the prominences in the experimental curve. 

The theoretical curves in figures 25 and 26 were calculated for 
T = 7800° K. Curve 1 involved a perturber density of N = N = 1.37 * 10"^ cm 
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Figure 2^.- Concluded. 
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folded with the Doppler profile for the 2 P - 6 D line. 
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Figure 25.- Concluded. 
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Figure 26.- Abel-inverted profile (r = 0) and quasistatic profiles folded 

3 3 

with the Doppler profile for the 2 P - 6 D line in the line core. 
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with no dependence of N on Av. Curve 2 involved N = 2 W with no depend- 
ence on Av. Curve 3 involved N = N with a dependence on Av of 

e 

1 



W) 


that is, a = 30 and b = 1. 

This choice of the parameters a and b improved the fit toward the 
center of the line in the violet wing, but no conclusions should be drawn here 
because of the uncertainties just mentioned. 

The Abel-inverted profile of the 2^P - 7 3 D line is presented with 
three theoretical profiles in figure 27* The parameters used in the theoret- 
ical calculations were: curve 1, a = 7, b = 1; and curve 2, a = 30, b = 1 

with N = N . Curve 3 has N = 2N^ and no dependence of N on Av was 
used. The fit in the wings is rather poor, mainly because of uncertainty in 
the background, although measurement inaccuracy is manifest. No improvement 

in the fit is indicated by using a wavelength-dependent perturber density. 

3 3 

In figure 28 the Abel-inverted profile of the 2 P - 8 D line is 

shown. In this case, the fit in the violet wing is quite good until well 

into the core of the line. The red wing fit is not good because of the 

presence of the forbidden component. Curve 1 was calculated with N = , 

a = 10, b = 2, curve 2 was calculated with N=N^, a = 10 , b=3, and 

curve 3 with N=N, a = 20, b = 2. Curve H was calculated with N = 2N 

e e 

and no dependence upon Av. It can be seen that the best overall agreement 
of the experimental profile was obtained with curve 3. It is possible, by 
further adjusting the parameters a and b, that even better agreement could 
be obtained. Several other values were tried, but the agreement was not as 


good. Those values are listed in table V. 
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TABLE V.- EMPIRICAL PARAMETERS USED TO EXTEND 
RANGE OF AGREEMENT OF QUASISTATIC THEORY 
FOR THE 2 3 P - 8 3 D TRANSITION 

N = 1.37 x 10 15 cm -1 
e 

T = 7800° K 


a 

b 

i 

1 

i 

2 

i 

10 

7 

1 

7 

10 

10 

1 

30 

1 

30 

2 

30 

10 

100 

1 

100 

2 

100 

10 


The following plotted in figure 28. 


7 2 

10 2 

10 3 

20 2 



xperime 



folded with the Doppler profile for the 2 P - 7 D line. 




Figure 27.- Concluded 





folded with the Doppler profile 
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In general, a fast transition from N = 2N^ to N = provided 
the best agreement. That this was so is indicated by having to go to larger 
values of b. The parameter a shifts the point of transition, whereas b 
governs the rate at which the transition occurs. 

The half widths of the averaged profiles and the Abel-inverted pro- 
files are compared with n(n - l) in figure 29* The points for the averaged 
profiles fall close to a straight line, as would be predicted by use of the 
hydrogenic quasistatic theory. The points for n = h and n = 5 fall above 
the line, indicating deviation from the quasistatic approximation and the 
influence of other types of broadening. The point for n = 10 has a large 
uncertainty because of overlapping lines and background. 

As was expected, the half widths corresponding to the Abel-inverted 
profiles (r = 0) lie above the points for the profiles averaged along the 
line of sight. These points also show the same trend with quantum number, 
except the point for n = 8 should be about 20 percent lower to fall on a 
line whose slope is the same as that of the averaged profiles. 

To determine whether one can use the half widths of the Abel-inverted 

2/3 

profiles to measure the electron density profile of the jet, ^i/2^ e 
compared with the radius in figure 30 . 

Note that the first few points of the n = 6 and n = 7 curves fall 
fairly well on a horizontal line; however, beyond that they begin to drop. 

This indicates that the latter two possibilities are more probable, because 
the reliability of the inversions decreases with increasing radius. The 
n = 8 curve begins to drop immediately and its relative position is too high, 
as indicated in figure 29. Thus, it is likely that in this case the half 



cm 


9h 



4 5 6 7 8 9 10 

n 

Figure 29.- Full widths of the measured Abel-inverted and non-Abel-inverted 

3 3 

line profiles of the 2 P - n D series. A best straight line is fit to the 
noninverted points. As expected the inverted data points lie above the 
noninverted points. The n = 8 point should fall about 20 percent lower 
in order for the Abel-inverted widths to have the same slope on the plot. 




0 0.2 0.4 0.6 0.8 1.0 

Normalized radius, r/R 


2/3 

Figure 30.- Full widths of the Abel-inverted profiles divided by are 

plotted against the radius. The decrease with radius indicates increasing 

2/3 

errors with r in the N or the widths or the failure of the N 

e 


dependence on width. 


9 6 


widths are in error as well, although it was in this case that the best agree- 
ment of the profile with the theoretical profile was obtained. Ideally, the 

curves should be straight horizontal lines. There are three possible causes 

2/3 

for the exhibited behavior: (l) the half widths do not obey the N 

e 

dependence; (2) the density profile is in error, with either the center values 
too low or the outer values too high; or (3) the half widths obtained from 
the Abel-inverted profiles are in error. To determine the probable cause of 
this behavior, the Abel-inversion technique must be investigated. It is 
possible that a better choice of expansion functions would improve the 
accuracy of the inversions. Probably, the Abel-inversions yield values too 
large toward the edge of the jet, thus causing the density to be too large. 



5 . 0 SUMMARY 


Until now, there have been virtually no helium line-profile measurements 

l4 15 -3 

m the electron density range from 10 to 10 cm ; all previous measurements 

13 3 16 IT 3 

have been made in low-density (10 cm or high-density (10 to 10 cm ~ 5 ) 

plasmas. This has prevented comparison of the quasistatic and impact theories 
for helium Stark broadening in a range where the uncertainties in the theories 
are greatest. The measurements presented herein indicate a general agreement 
with the quasistatic theory for electrons and ions in the wings, except where 
forbidden components or overlapping with neighboring lines prevent a valid 
comparison. There was a general improvement in the agreement of the calculated 
profiles with the measured ones when the dependence of perturber density on 
wave number was included empirically in the calculations. This is consistent 
with the results of Schliiter (ref. 7 ), and an even better consistency with 
his results is found if Unsold's estimate (ref. 26) for the critical value for 
the transition between the quasistatic and impact domains of electron broadening 
Av^ is used. His estimate is approximately one-half that of Traving (ref. 19); 
using it would effectively decrease the parameter a by a factor of one-half, 
yielding values closer to those of Schliiter in his hydrogen-profile calcula- 
tions . 

Because the only method used for determining the electron density was 
absolute continuum measurements , no assessment of any systematic errors could 
be obtained; and, because the electron density is the most important parameter 
in the profile calculations , future comparisons should have the benefit of a 
more certain electron density. 
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To improve the accuracy of the measurements, better accounting of the 

-5/2 

background is needed. If the assumption is made that the wings obey the Av 

dependence , the fit in the wings could be forced by subtracting an amount of 

-5/2 

background which will yield the best Av fit. Also, a more precise Abel- 

inversion measurement and numerical technique should be attempted. Choosing 
expansion functions which are a set of orthogonal functions satisfying suitable 
boundary conditions, as suggested by Birkeland and Oss (ref. 27), might improve 
the numerical fit in the Abel-inversion. 

Because the standard error in is ±13 percent at the jet center 

line and the line profile measurements have about ±5 to 10 percent uncertainty 
(possibly larger in the far wings), the comparisons have an accuracy of ±13 to 
l6 percent, at best. 

To achieve the best overall agreement in the line core and in the wings , 
the impact theory would have to be accounted for in the theoretical profiles. 
However, this was not the objective of this work. This may be done in future 


investigations . 


APPENDIX A 


ABSOLUTE CALIBRATION OF THE SPECTROMETER 

The response of a voltmeter which measures the current in the spec- 
trometer photomultiplier is the folding of the instrument profile with the 
source profile. This can he represented mathematically by the equation 


J .oo 

- O 1 


(X) dA 


(Al) 


where V is the measured voltage, l(X) is the source intensity profile, 

a is a normalization constant, and P(X) is the normalized shape function 
A 

r 

of the instrument profile, | P(X) dX = 1. 

•'0 


Suppose that the source is a calibrated source, such as a ribbon- 
filament lamp, whose intensity function is denoted by I (X). If it is 

CEl-L 

assumed that I (x) is very nearly constant in an interval in which P(X) 
c a_L 

is not negligible (X < X < X ^) , then equation (Al) will reduce to 



(A2) 


where I , « I ..(X ) and V . is the voltmeter response to the calibra- 
cal cal o cal 

tion source. Because of the normalization condition, equation (A2) will re- 
duce to V , = I , o. . If the sensitivity of the spectrometer S, is 
cal cal X X 

defined by the relation V = SI, where V is the response of the voltmeter 

A 
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to an intensity I, then it is readily seen that the normalization constant 

a. equals the sensitivity S . The absolute sensitivity is then 
A A 


a 


X 




cal 

cal 


where I is in absolute units, 
cal 


(A3) 


Continuum Radiation 

To measure the steradiancy of an unknown continuum source, the source 
is placed in the same location as the calibration source and the voltage V c 
is measured 


V 

c 


V 


cal^c 

■'"cal 


(aU) 


Then, intensity of the unknown continuum source in absolute units is 


I 

c 


I cal V c 


V 


cal 


V 

c 



(A5) 


Line Radiation 

If the source is line radiation rather than continuum, the treatment 
is somewhat more involved. There are two cases. One case is where the width 
of the instrument profile is not wide compared with the line to be measured. 

In this case, one must scan over the line and integrate the resulting profile. 
The second case is where the instrument profile is wide compared with the 
line. If background is negligible, then only the peak (flat top) intensity 


must be measured. 
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Case I. Instrument profile not -wide compared with the line profile . - 
The output of the voltmeter can he represented by the equation 


v( A ' ) 



cr, P( A - A ' )<J ( A ) dA 
A 


(A6) 


where J(A) is the spectral steradiancy of the line in yW/cm -sr-A. 
If 


I 


line 



j( A ) dA 


(at) 


is the total steradiancy of the line in yW/cm -sr, 

I. . L(A) = J(A) and I L(A) dA = 1, where L(A) 
line J 


profile . 


then we can write 
is the normalized line 


By scanning over wavelength then integrating the result, we have 


V 


line 




)l_, . L( A )dAdA ' 
line 


(A8) 


O 

V, . has units of mV - A, say. 
line 

Because P(A) and L(A) are independent, the order of integration 
can be interchanged. 



A ' ) dA ' = I 


line 



P( A - A ' )dA ' 
A 



L(A)dA 


(A9) 


Then , 


V 


line 


I, • s 
line 


A* 
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Solving equation (A9) for I an< ^ substituting equation [Ah) 

Cor S x , 


line 


1 1 V T 

cal line 
V cal 


(A10) 


in uW/sr-cm . 

Case II. Instrument profile wide compared with the line profile .- 
Beginning with equation (A6), the assumption is being made that P(a) is 
slowly varying or is constant over the range where j(A) is nonvanishing. 
Then, it can be removed from the integral 


v( A ' ) 




J(A) dA 


(All) 


where P^° is the peak value of the normalized profile. By equation (A7 ) , 
this becomes 



P, a I . 

A A line 


(A12) 


Using equation (A6) and integrating (scan) over the profile. 


Jv(X')dX' = ffc 


(A')dA' = lo. P(A - A ' ) J( A )dAdA 1 


= I, . la. P(A - A' ) dA 
line J A 


(A13) 


is obtained. 
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Dividing equation (A12) by equation (A13) 




j"v( A ' )dX 


p a° ct a 

cj^P( A - A ' )dA 


(AlH) 


Using equation (A3) and solving for 


°A P A 


J 




v( A 1 )dA 1 


(A15) 


is obtained. 

Knowing A o P^° f or an y line, equation (A12) can be used to calibrate 
any narrow line by a single measurement at the peak intensity v . Therefore, 


v 

I = fe- 
line 




v( A )dA 

IkT 


(A16) 


In the presence of background, the background voltage must be sub- 
tracted before integrating (if integration over A is required). Otherwise, 
v^(A^ - A^) can be subtracted from the integrated value, where A^ and A^ 
are the lower and upper limits of the integration, and v^ is the background 
signal . 

The units of Jv( A)dA/v(A Q ) are wavelength. This quantity can be 
defined as the effective spectral slit width W. If the profile were an 
ideal rectangle, triangle, or trapezoid, then W would be the full half-width 
of the profile. When the exit slit is wide compared with the entrance slit, 

W is very nearly equal to the spectral exit-slit width. 


APPENDIX B 


ABEL-INVERSION 

The Abel-inversion technique that was used in these experiments was 
adapted from that of Freeman and Katz (ref. 25). A review of the technique is 
given here. 

Consider a cylindrically symmetric luminous plasma, viewed perpendicu- 
larly to the flow axis. The plasma is assumed to he optically thin. The ob- 
served intensity is the sum of intensities emitted from volume elements along 
the line of sight (fig. B-l). 

These measurements are essentially averages or integrated intensities 
over the thickness of the plasma, which is assumed to have cylindrical symmetry. 
If the contribution to the intensity at a radius r is N(r), the measured in- 
tensity can be represented by 



(Bl) 


(fig. B-l). The integration is along the line of sight. 


Changing the integration variable to r yields 



(B2) 


where x 


( 2 2 

\r - y 



10k 
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Equation (B2) is a form of Abel's integral equation, and it has the 
formal solution 


N(r) 


lid 
* r dr 




yi(y) 



dy 


(B3) 


If l(R) = 0, then equation ( B3 ) can he written by integration by 

parts as 


N(r) 


1 f K dl(y) ( 2 2\ 1/2 , 

s -*l ~ r > dy 


(B4) 


Unfortunately, the analytical form of l(y) is not known, so equation (b 4) 
cannot be integrated directly. A numerical technique must be used. One method 
is to expand h(r) or l(y) in a set of functional relations, then carry out 
the integration of this expansion analytically. Thus, if N(r) is approxi- 
mated by a linear combination of a set of functions 


N(r) 




r) 


(B5) 


then l(y) =7 a n, (y) where E, (r) and the Abel transforms of n, (y) . 

jl k k k k 

The expansion coefficients a are determined by fitting experimen- 
tal data J(y ), which is the measured intensity as a function of position, 
to a set of expansion functions by the method of least squares. Then, having 
these expansion coefficients and the analytically transformed expansion func- 
tions £ (r) , N(r) can be calculated easily. 

K 
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The assumption of perfect cylindrical symmetry of the experiment is 
not entirely justified. It then becomes convenient to symmetrize the data 
as follows . 

J + (y) = 77 [J(y) + J(-y)] symmetrical part (b6) 

J (y) = 7 T [J(y) - j(-y)] antisymmetrical part (B 7 ) 


In the case of the Abel-inverted line profiles, each of these parts 
was fit by a least-squares technique to functions of the type 


I+(y) - Ev/w 


(B8) 


r<y) = E\"k' (y) 


(B9) 


-J. _ 

where (y) are symmetric functions and n k ~(y) are antisymmetric func- 
tions . 

A measure of the goodness of fit is given by the error formula 


a 




(BIO) 


The number of terms in the expansion was chosen such that a was a 
minimum. The Abel-inverted functions are given by 


N+(r) = ZA ? k +(r) 


(Bll) 
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N (r) = T^ Vk (r) (B12) 

where £ + (r) and £ (r) are the transforms of n, + (y) and n, (y) , respec- 

K. K K K 

tively . 

The following functions were used. 


n k +( y) = (r 2 - y 2 ) 




(B13) 


s k + (r) = x k (e 2 
S k -(r) = */(** 


k-1/2 


k-1/2 


where 


X 


k 


1 p2k (k!) 2 

tt (2k) ! 


(BlU) 


(B15) 


( Bl6 ) 


A computer program was written to compute the expansion coefficients 
a k and t> k By the least-squares method. Then, the program was used to com- 
pute N (r) compared with r for a set of radii. 

In the cases other than the Abel-inverted line profiles, the symmetri- 
zation just described was not done. Here the data J(y^) were fitted to the 
symmetric functions, and no antisymmetric part was calculated. 
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